从分子水平上研究群体的遗传结构特征，不仅可以了解群体遗传分化水平，同时也有利于遗传育种和优良群体遗传保护。线粒体DNA(Mitochondrial DNA,  mtDNA)与核DNA相比，碱基排列与堆叠方式相对简单、反映母系遗传特征，传代过程中变异速度较快、核苷酸发生转换或者颠换的频率较高等，mtDNA的这些独特之处使其成为鱼类群体遗传学研究中被众多研究人员广泛应用的重要标记方法(Oleinik et al., 2007; Teletchea, 2009)。mtDNA基因组内不同区域的碱基替换呈现出不同规律和特征，适合应用于不同目标和层次群体遗传研究。在mtDNA上，控制区(D-loop)是整个mtDNA上序列和长度变异最大的区域(Lee and Kocher, 1995)，也是整个mtDNA基因组中进化最快的区域，通常适用于种内群体间遗传差异检测和分析，也有用于种间遗传差异分析的研究报道(郭新红等, 2004; 颉晓勇等, 2011)。细胞色素Cytb基因是蛋白编码基因，其进化速度与mtDNA其它区域相比较处于中间水平，适合种群间遗传差异的检测，是在DNA分子水平上分析生物种间和种内遗传变异的良好方法(Teletchea, 2009)。然而，有关Cytb和D-loop序列遗传变异信息研究结果之间区别和联系的研究报道相对较少，特别是对于两种技术所得到的遗传变异数据之间关系所知无几。

吉富罗非鱼由上海水产大学引进中国之后，经历10余年时间的遗传改良，成为生产性状优良的养殖新品系(颉晓勇等, 2011)，对选育群体从分子水平上进行深入的遗传分析具有重要的意义。本研究开展了新吉富罗非鱼选育群体mtDNA Cytb和D-loop基因序列变异信息的对比分析，旨在探索2种mtDNA标记分析结果之间关系，目的是为了针对不同研究材料，选择不同的分析方法，以及扩大研究结果之间横向比较在分析技术上提供借鉴和参考，同时使得相关研究分析更加全面、研究结论更加客观。

1  结果与分析

1.1  基于Cytb与D-loop序列分析的群体内遗传多态性比率

本研究基于Cytb与D-loop序列分析得到的群体内遗传多样性比率结果如表1所示。根据2种方法得到的多态位点比率范围为0.67905~0.80492，平均0.77674；单倍型多样度比率处于0.86412~1.07612之间，平均0.91945；核苷酸多样性指数比率处于0.71377~0.88449之间，平均0.76977；平均核苷酸差异数比率处于0.89954~1.08212之间，均值为0.93619。

1.2  Cytb与D-loop序列中碱基转换率和颠换率

本研究Cytb与D-loop序列中碱基转换率和颠换率对比如表2所示。Cytb序列碱基转换率处于0.03195~0.05921之间，平均0.04267；Cytb碱基颠换率处于0.00282~0.00548之间，平均0.00410。D-loop序列碱基转换率处于0.02772~0.04767之间，平均0.03725；D-loop碱基颠换率处于0.01552~0.02772之间，平均0.02284。

2  讨论

在遗传方式方面，mtDNA以非孟德尔遗传方式控制和影响线粒体功能，与核基因不相一致。同时，鱼类mtDNA的碱基突变率高于核内DNA，mtDNA损伤缺少有效的修复机制，特别需要指出的是有些mtDNA突变还存在累加效应，对mtDNA多态现象的量化分析表明两个无关个体之间mtDNA碱基平均相差3%(郭新红等, 2004)。因此，对群体进行严谨的遗传变异特征分析时，有必要选择合适的mtDNA区域进行分析。有大量采用Cytb、D-loop或其它mtDNA基因区域进行分析的研究报道(陈善元和肖蘅, 2003; 于旭蓉等, 2011; Gerlach and Musolf, 2000)。mtDNA不同区域的进化速率不同，的研究报采用mtDNA不同区域分析得到的序列遗传变异特征必然有所差异。但是目前对于mtDNA不同区域遗传得到的序列变异特征进行比较和量化分析道较少，当前分子标记相关研究多数停留于单种标记分析而缺乏不同种类标记之间横向分析，较大程度上限制了研究结果的参考价值，因而这方面仍有大量工作有待深入研究。

mtDNA序列碱基变异可以分为两种类型，转换(transition)是发生同种类碱基置换，即由嘌呤或者嘧啶发生同种类但是不同碱基之间的置换；颠换(transversion)是发生不同种类碱基之间的置换，嘌呤置换嘧啶或嘧啶置换嘌呤，即发生嘌呤和嘧啶碱基之间的置换。本研究Cytb转换率平均0.04267，Cytb颠换率平均0.00410；D-loop转换率平均0.03725，D-loop颠换率平均0.02284。可以得到Cytb区域转换/颠换比值约为10.41，D-loop区域转换/颠换比值约为1.66，转换率均大于颠换率，该结论与与其他的研究结果相一致(Gerlach and Musolf, 2000; 杨慧荣等, 2012; 赵亮等, 2010)。碱基置换在不同基因区域具有不同生物学功能，发生于编码多肽的基因区域比如Cytb，则可能改变密码子信息从而使得后续转录、翻译等相应改变，生物体中可能出现一种新的氨基酸结构取代相应位置原有的某一种氨基酸；碱基置换也可能使蛋白编码基因中提前出现终止密码，多肽链合成提前中断。这些碱基置换的结果都使得蛋白编码基因不能正确地形成原有的蛋白质，从而完全或者部分地失去某种生物学活性(贾若欣等, 2011, 畜牧与兽医, 43(9): 100-104)。本研究不同基因区域碱基转换率Cytb/D-loop约为1.15，不同基因区域碱基颠换率Cytb/D-loop约为0.18，推测该结果与Cytb编码蛋白，而D-loop与mtDNA转录和复制过程中的调控作用有关。

DNA双链中的不同区域具有不同突变频率，以及承受不同的进化压力，因而具有不同的进化速率，本研究也证实了这一结论。D-loop区域的变异速率约为mtDNA完整分子的5~10倍(郑冰蓉等, 2002)。本研究尼罗罗非鱼5个选育群体Cytb/D-loop序列多态性比率结果表明，多态位点比率和核苷酸多样性指数Cytb仅仅是D-loop区域的0.77674和0.76977，Cytb/D-loop单倍型多样度和平均核苷酸差异数比率略高，分别为0.91945和0.93619，这一结果与控制区是mtDNA中不编码多肽链的核苷酸片段，无修复系统，不受选择压力的影响，因而积累了更多的变异(李鹏飞等, 2011)，而细胞色素b基因编码蛋白因而进化相对较慢有关。总体来说，本研究对Cytb和D-loop区域序列变异信息的对比分析，有助于为其它类似研究扩大横向比较范围和更深入地理解单种Cytb或D-loop标记分析方法的研究结果提供参考。

3  实验材料与方法

3.1  实验材料

上海水产大学1994年从菲律宾引进的5000尾尼罗罗非鱼(基础群体, F₀)，从1996年起开展选育，每年产生一代。本实验的材料是从F₀、F₆－F₉共5代尼罗罗非鱼选育群体中随机采样各20尾，剪取尾鳍，分别编号后以95%乙醇保存备用。

3.2  基因组DNA提取和PCR扩增

参照常规的酚－氯仿抽提程序进行，并通过琼脂糖凝胶电泳初步检测所提取DNA质量。控制区引物序列参考Agnèse等(1997)，碱基序列为 DL1：5’ ACC CCT GGC TCC CAA AGC 3’ 和 DH2：5’ ATC TTA GCA TCT TCA GTG 3’。细胞色素b基因引物参考彭作刚等(2005)，碱基序列为 L14724：5’ GAC TTG AAA AAC CAC CGT TG 3’ 和 H15915：5’ CTC TCT CCG GAT TAC AAG AC 3’。PCR扩增反应体系为50µL，包含10×Buffer 5µL，dNTP 10 mmol，上游引物和下游引物各2 µL (10 pmol/µL)，模板DNA 200 ng，Taq DNA聚合酶2.5U。PCR扩增程序：94℃预变性设定4min；94℃高温变性设定30s，模板与引物适温条件下退火30s，72℃引物聚合与延伸1 min，每次扩增程序包含35个循环；最后72℃10 min；4℃保存。其中，控制区退火温度50℃，细胞色素b基因退火温度56℃。

3.3  DNA测序

扩增产物用1%的琼脂糖凝胶电泳进行检测，选择扩增效果良好的PCR产物进行纯化，由上海sangon公司用ABI 377 DNA自动测序仪进行双向序列测定。

3.4  数据处理与分析

用BLAST软件(Altschul et al., 1997)搜索GenBank中尼罗罗非鱼的控制区和细胞色素b基因序列，与实验测定的尼罗罗非鱼相应序列结果进行比较；用Bioedit软件(Hall, 1999)编辑基因测序结果，以CLUSTL W软件(Thompson et al., 1994)重排序列和比较同源性，并辅以人工核对校正。用DNASP软件(Rozas et al, 2003)统计分析。按Cytb/D-loop分别计算由Cytb和D-loop 2种方法所得到的多态位点比率、单倍型多样度比率、核苷酸多样性指数比率、平均核苷酸差异数比率。按转换与颠换数占相应序列碱基总数比率计算转换率与颠换率。

参考文献 

Agnèse J.F., Adépo-Gourène B., Abban E.K., and Fermon Y., 1997, Genetic differentiation among natural populations of the Nile tilapia Oreaochronis niloticus (Teleostei Cichlidae), Heredity, 79: 88-96

http://dx.doi.org/10.1038/hdy.1997.126

Chen S.Y., and Xiao H., 2003, Molecular markers and inference methods in molecular phylogenetics of fishes, Yunnan Daxue Xuebao (Ziran Kexue Ban) (Journal of Yunnan University (Natural Sciences)), 25(Suppl.): 146-152 (陈善元, 肖蘅, 鱼类分子系统发育研究中的分子标记与推断方法, 云南大学学报(自然科学版), 25(增刊): 146-152)

Gerlach G., Musolf K.F., 2000, Fragmentation of landscape as a cause for genetic subdivision in bank voles, Conservation Biology, 14(4): 1066-1074

http://dx.doi.org/10.1046/j.1523-1739.2000.98519.x

Guo X.H., Liu S.J., Liu Q., and Liu Y., 2004, New progress on mitochondrial DNA in fish, Yichuan Xuebao (Acta Genetica Sinica), 31(9): 983-1000 (郭新红, 刘少军, 刘巧, 刘筠, 2004, 鱼类线粒体DNA研究新进展, 遗传学报, 31(9): 983-1000)

Hall T.A., 1999, BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT, Nucleic Acids Symposium Series, 41: 95-98

Lee W.J., Kocher T.D., 1995, Complete sequence of a sea Lamprey (Petromyzoa marinus) mitochondrial genome: Early establishment of the vertebrate genome organization, Genetics, 139(2): 873-887

Li P.F., He Z.T., Xu K.D., Zhu W.B., Xue L.J., and Zhang H.L., 2011, Genetic diversity of Pampus cinereus in the east China sea revealed by mtDNA D-loop sequence, Zhejiang Haiyang Xueyuan Xuebao (Journal of Zhejiang Ocean University), 30(1): 14-17 (李鹏飞, 贺舟挺, 徐开达, 朱文斌, 薛利建, 张洪亮, 2011, 东海野生灰鲳mtDNA D-loop序列遗传变异分析, 浙江海洋学院学报(自然科学版), 30(1): 14-17)

Oleinik A.G., Skurikhina L.A., and Brykov V.A., 2007, Divergence of Salvelinus species from northeast Asia based on mitochondrial DNA, Ecology of freshwater Fish, 16(1): 87-98

Peng Z.G., Zhang Y.G., He S.P., and Chen Y.Y., 2005, Phylogeny of Chinese catfishes inferred from mitochondrial cytochrome b sequences, Yichuan Xuebao (Acta Genetica Sinica), 32(2): 145-154 (彭作刚, 张耀光, 何舜平, 陈宜瑜, 2005, 从细胞色素b基因序列变异分析中国鲇形目鱼类的系统发育, 遗传学报, 32(2): 145-154)

Rozas J., Sánchez-DelBarrio J.C., Messeguer X., and Rozas R., 2003, DnaSP, DNA polymorphism analyses by the coalescent and other methods, Bioinformatics, 19(18): 2496-2497

http://dx.doi.org/10.1093/bioinformatics/btg359

Altschul S.F., Madden T.L., Schäffer A.A., Zhang J., Zhang Z., Miller W., and Lipman D.J., 1997, Gapped BLAST and PSI-BLAST: A new generation of protein database search programs, Nucl. Acids. Res., 25(17): 3389-3402

Teletchea F., 2009, Molecular identification methods of fish species: Reassessment and possible applications, Reviews in Fish Biology and Fisheries, 19(3): 265-293

Thompson J.D., Higgins D.G., and Gibson T.J., 1994, CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through weighting, positioin-specific gap penalties and weight matrix choice, Nucl. Acids. Res., 22(22): 4673-4680

http://dx.doi.org/10.1093/nar/22.22.4673

Xie X.Y., Li S.F., and Cai W.Q., 2011, Analysis of genetic diversity of tilapia during selection processing based on D-loop sequence, Shanghai Haiyang Daxue Xuebao (Journal of shanghai ocean university), 20(3): 336-341 (颉晓勇, 李思发, 蔡完其, 2011, 基于D-loop序列的罗非鱼选育群体遗传变异分析, 上海海洋大学学报, 20(3): 336-341)

Yang H.R., Zhao H.H., Mong Z.N., Liu L., and Lin Q.Z., 2012, Comparison of mitochondrial D-loop and cytb sequences of Squaliobarbus curriculus, Zhongshan Daxue Xuebao (Ziran Kexue Ban) (Acta Scientiarum Naturalium Universitatis Sunyatseni), 51(5): 100-106, 136 (杨慧荣, 赵会宏, 蒙子宁, 刘丽, 林权卓, 2012，赤眼鳟线粒体D-loop和Cytb基因序列的对比分析, 中山大学学报(自然科学版), 51(5): 100-106, 136)

Yu X.R., Qiu X.M., Liu X.Y., and Xu L., 2011, Application of research on marine animals population genetic structure using mitochondrial DNA polymorphism, Shengwu Jishu Tongbao (Biotechnology Bulletin), 10: 49-54 (于旭蓉, 仇雪梅, 柳晓瑜, 徐丽, 2011, 线粒体DNA多态性在海洋动物群体遗传结构研究中的应用, 生物技术通报,  10: 49-54)

Zhao L., Xie B.G., Liu Z.J., Xu M.Q., and Li M., 2010, Molecular structure and DNA substitution rate of the mitochondrial control region and cytochrome b in Taihu Salangid, Neosalanx taihuensis, Dongwuxue Zazhi (Chinese Journal of Zoology), 45(2): 27-38 (赵亮, 谢本贵, 刘志瑾, 许木启, 李明, 2010, 太湖新银鱼线粒体D-loop和cytb片段序列结构与进化速率比较, 动物学杂志, 45(2): 27-38)

Zheng B.R., Zhang Y.P., Xiao H., Lan J.H., and Zan R.G., 2002, The sequence variation feature of mt DNA D-loop region of Cyprinus, Shuichan Xuebao (Journal of fisheries of China), 26(4): 289-294 (郑冰蓉, 张亚平, 肖蘅, 蓝家湖, 昝瑞光, 2002, 鲤属鱼类mtDNA控制区(D-环区)序列的变异性分析, 水产学报, 26(4): 289-294)